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Abstract 



The partial decay rate of the Z boson into bottom quarks constitutes an important 
decay channel. This is mainly due to the virtual presence of the top quark in the 
loop diagrams giving rise to correction factors which are quadratic in the top quark 
mass. At one- and two-loop order it turned out that the leading term in the heavy- 
top expansion leads to very good approximations to the exact result. In this work 
w ■ the non-singlet diagrams at 0{o? s GfM^) are considered. 



The impressive experimental precision mainly at the Large Electron Positron collider 
(LEP) at CERN, the Stanford Linear Collider (SLC) and the FERMILAB Tevatron in 
Chicago has made it mandatory to evaluate higher order quantum corrections to the 
processes observed in the experiments |I[. The strategy to combine experimental infor- 
mation with theoretical computations has successfully been applied to the search for the 
top quark several years ago. Nowadays the same concept is used in order to pin down the 
mass of the Higgs boson, the only not yet discovered particle of the Standard Model of 
elementary particle physics. 

An important observable is the decay of the Z boson into bottom. QCD corrections 
are known up to 0(a^) (for a comprehensive review see @). The electroweak one-loop 
corrections are known since quite some time ||. They have the interesting feature that 
the top quark appears virtually in the loop diagrams. Recently also the full corrections of 
0(aa s ) were completed @, |], ||. The diagrams involving a top quark are considered in |J 
where the first five terms in the expansion for a heavy top quark mass, M t , is computed. 
It was demonstrated that these terms approximate the exact result quite well. Actually 
it turned out that both at 0{a) and 0(aa s ) a large cancellation between the sub-leading 
terms takes place and effectively only the leading term proportional to GfM^ |7j remains. 
This is a strong motivation to look at the next order in the strong coupling constant and 
evaluate the leading terms. To the corrections enhanced by the top quark mass only those 
contributions have to be considered where a scalar particle, namely the Higgs boson, H, 
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the neutral Goldstone boson, x, or the charged one, </>*, couples to the top quark. Thus 
no diagrams have to be considered where the W or Z boson appear as internal lines. 

Corrections of this order were first computed for the p parameter ||, the ratio of the 
charged and neutral current amplitude, where it turned out that they are quite impor- 
tant |§. Later on also the hadronic Higgs decay was analyzed at (9(a;?GVM 4 2 ) |I| |Tl| . 
In the case of the Higgs boson one can exploit that the scalar coupling is proportional to 
the mass which simplifies the construction of an effective Lagrangian and especially the 
subsequent evaluation of the diagrams. Actually the whole computation could be reduced 
to the evaluation of two-point functions. We will see below that in the case of the Z boson 
one should also consider vertex diagrams. 

It has become customary to parametrize the corrections proportional to M 4 2 by the quatity 
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respectively the quantity x t which is defined using the MS definition of the top quark 
mass, m t . 

The quantum corrections to T(Z — > bb) are divided into universal ones which are identical 
for all quark species and non-universal parts which are specific for the Zbb vertex. Both the 
universal and non-universal corrections are divided into singlet and non-singlet parts. The 
singlet contributions arise from those diagrams where the Z boson and the bottom quarks 
of the final state couple to different fermion lines. Another class of singlet contributions is 
constituted by the diagrams where the Z boson couples to two charged Goldstone bosons 
which in turn form together with two gluons a box diagram and the gluons finally couple to 
the quarks in the final state. In Fig. [I] some sample diagrams are listed. Fig. [l|(a) and (b) 
are of universal nature whereas in (c) the Goldstone boson in directly coupled to the final 
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Figure 1: Singlet diagrams contributing at 0{a 2 s X t ) to the hadronic Z boson decay. In 
(a) and (c) the dashed line correspond to the charged Goldstone boson whereas in (b) 
it may also be the Higgs or the neutral Goldstone boson. Diagrams (a) and (b) are of 
universal type whereas (c) constitutes a non- universal contribution to T(Z — > bb). In 
the displayed examples the thick lines correspond to top quarks whereas the thin lines 
represent bottom quarks. 
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Figure 2: Non-singlet contributions of (9(a^X 4 ) to the hadronic Z boson decay. In (a) 
the dashed line corresponds to the Higgs boson or the neutral or charged Goldstone boson. 
In the diagrams (b) — (e) only the charged Goldstone boson is allowed. Diagram (a) is 
of universal type whereas (b) — (e) constitute non-universal contributions to T(Z — > bb). 
In the displayed examples the thick lines correspond to top quarks whereas the thin lines 
represent bottom quarks. 



state bottom quark thus providing a non-universal contribution. In this article only non- 
singlet diagrams will be computed; the singlet contributions will be considered elsewhere. 
The universal corrections of 0(a 2 s X t ) are in part governed by the p parameter ||. A second 
source for universal corrections arise from those diagrams where only gluons couple to the 
light quark lines. The gluons split into a fermion loop actually formed by bottom and 
top quarks accompanied by an additional exchange of a scalar particle (cf. Fig. 0(a)). 
The main focus of this paper is devoted to the evaluation of the non-universal non-singlet 
diagrams. Typical examples are pictured in Fig. 0(o)- (e). 

In a first step an effective Lagrangian is constructed were the top quark is integrated 
out. Thereby it is convenient to split the fermion fields into their left and right part and 
consider them separately. It is furthermore necessary to decouple the bottom quark fields 
using the relations (see e.g. 
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where the primes denote the quantities in the effective theory and the superscript "0" 
reminds that we are still dealing with bare quantities. The decoupling constants ( 2 'b 



can be computed with the help of 
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where E" y and £" a is the vector and axial- vector part of the bottom quark self energy. 
Here only the hard part, i.e. those diagrams containing the top quark, has to be computed 
which is indicated by the index "h" . Finally the part of the effective Lagrangian describing 
the interaction of the Z boson to bottom quarks has the form 
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The residual dependence on M t is contained in the coefficient functions C[, R . They are 
obtained from the hard part of the Zbb vertex: 
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Here, the left and right part of the Zbb vertex are defined through: 
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where it is understood that in addition to the top-induced diagrams also the tree-level 
terms are included. The coefficient functions of Eq. (|5|) are finite after the coupling con- 
stant a s and the mass of the top quark are expressed through their renormalized coun- 
terparts. This is because the vector and axial-vector currents have vanishing anomalous 
dimension as long as only non-singlet diagrams are considered. Thus from now on the 
index "0" is omitted. 

In order to evaluate the partial decay rate of the Z boson into bottom quarks at 0(a 2 s X t ) 
one has to evaluate the coefficient functions up to this accuracy. Furthermore the pure 
QCD corrections in the effective theory are needed up to order a 2 s . It can be taken over 
from [O] and reads: 
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where rii 
value ( 3 p 



= 5 is the number of light quarks and Q is Riemann's Zeta function with the 
1.202056903. 



The computation of the decoupling constants for the bottom quark field up to order 
<y 2 X t has been performed in JL(|. The only missing pieces are the vector and axial- vector 



contributions to the hard part of the Zbb vertex. Some sample diagrams are listed in 
Fig. 0. As mentioned above only those diagrams have to be taken into account which 
contain a virtual top quark. Note that for the very calculation it is possible to nullify 



all external momenta. At one-loop order only two diagrams have to be considered. This 
increases to 19 at the two-loop level which is still feasible by hand. In the order we 
are interested in, however, more than 350 diagrams have to be considered, which makes 
the extensive use of computer algebra necessary. For the present calculation the package 
GEFICOM [[0J has been used. It passes the generation of the diagrams to QGRAF |TJ]] and 
uses for the very computation of the integrals the program MATAD |15j which is written 
in FORM |16|] for the purpose to compute one-, two- and three-loop vacuum graphs. For a 



recent review concerned with the automatic computation of Feynman diagrams see [[17 
Expressed in terms of the MS top quark mass the result for the coefficient functions read: 
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with m t = m t (fi). After the second equal sign the colour factors Cp = 4/3, C 'a = 3 and 
T = 1/2 have been inserted. (2 = vr 2 /6 and £4 = 7r 4 /90. S0 and ce are the sine and cosine 
of the weak mixing angle. The constants B4, D3 and S2 typically appear in the result of 



three- loop vacuum integrals and read [18 
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Note that according to the QED Ward identity the universal corrections induced by the 
diagrams in Fig. [|(a) cancel in the coefficient functions against the corresponding part in 
the quark self energy. As we consider in addition the bottom quark to be massless the 
right-handed coefficient function sticks to its Born value. Using the relation between m t 
and the on-shell mass M t [2(J one gets: 
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Let us now turn to a brief numerical discussion of the new results. The decay rate can be 
computed with the help of 
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Actually two scales are involved in the process, namely Mz and the mass of the top 
quark. The resummation of potentially large logarithms is, however, trivial as both Cl/r 
and S^'^ CD are separately renormalization group invariant. Thus the scale parameter /x 
may be set to m t , respectively, M t in the coefficient functions and to Mz in the massless 
corrections. For these choices the numerical expansions of the ingredients for Eq. flTT|) 
read: 
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with /i t 

Concerning the enhanced corrections of 0(X t ) to the coefficient functions the same ob- 
servations can be made as for the p parameter || and the various quantities in connection 
with the Higgs decay [^1, 10 1: Expressed in terms of the on-shell top quark mass the 
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leading order term is "screened" by the QCD corrections as they enter with a different 
sign. On the other hand, the coefficients turn out to be much smaller in the MS scheme. 
Actually the coefficient in front of the three-loop term is smaller by a factor of 50 as com- 
pared to the corresponding one in the on-shell scheme. Furthermore the sign is alternating 
which also indicates a faster convergence if the MS mass is used for the parameterization. 

Inserting Eqs. (|T^) into Eq. ( |TT| ) finally leads to the following M^-enhanced terms: 

r**(Z->&6) = N f%\ a (-1 + \si) x t [1 + 0-0161 + 0-0014] , (13) 

r*(Z-65) = ^^(-1 + LgW [1-0.074 -0.0130], (14) 

where the values af\M z ) = 0.118, af\M t ) = 0.107 and af\m t ) = 0.108 have been 
used. The numbers in the squared brackets correspond to the corrections with increasing 
power in a s . The index OS reminds on the definition of the top quark mass used and as 
before the index MS is suppressed. In the MS scheme the second order QCD corrections 
amount to roughly 9% of the first order ones, however, the overall size is quite small. The 
order a s corrections to the leading X t term in the on-shell scheme is almost by a factor 
of five larger than in the MS scheme and the corresponding 0(a 2 ) corrections amount 
to almost 18% of the 0(a s ) term. It is actually almost as large as the 0(a s ) term in 
Eq. ©. 

To summarize, quantum corrections of 0(a 2 Xt) to the decay of the Z boson into bottom 
quarks have been computed. If we assume that the observations made at one- and two- 
loop level are also true at order a 2 X t a substantial part of the corrections is available. 
Expressed in terms of the MS mass they turn out to be tiny. In the on-shell scheme the 
quantum corrections are much larger and they screen the leading X t term by almost 9%. 
Note that the newly computed term of 0(a 2 x t ) makes it possible to use the combination 
of the three- loop p parameter || and the partial width T(Z — > bb) in a consistent way. 
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